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The two-parameter Poisson-Dirichlet distribution is a probabihty distribution on the totality of 
positive decreasing sequences with sum 1 and hence considered to govern masses of a random 
discrete distribution. A characterization of the associated point process (that is, the random 
point process obtained by regarding the masses as points in the positive real line) is given in 
terms of the correlation functions. Using this, we apply the theory of point processes to reveal the 
mathematical structure of the two-parameter Poisson-Dirichlet distribution. Also, developing 
the Laplace transform approach due to Pitman and Yor, we are able to extend several results 
previously known for the one-parameter case. The Markov-Krein identity for the generalized 
Dirichlet process is discussed from the point of view of functional analysis based on the two- 
parameter Poisson-Dirichlet distribution. 

Keywords: correlation function; Markov-Krein identity; point process; Poisson-Dirichlet 
distribution 

1. Introduction 

This paper is concerned with a two-parameter family of probability distributions on the 
infinite-dimensional simplex Voo of non- negative decreasing sequences with sum 1, 



\{vi) = {vi,V2, . . .) > z;2 > ■ ■ ■ > 0,y^«i =^ l|- 



It extends the one-parameter family of distributions known as Poisson-Dirichlet dis- 
tributions, which were introduced by Kingman [32]; see, for example, [1, 44, 50] for 
related topics and bibliographic information. Pitman and Yor [45] defined the two- 
parameter Poisson-Dirichlet distribution, denoted PD(a,0), in the following manner. 
Given < a < 1 and 6 > ~a, define a sequence (Vi) of random variables by 

Vi^Yi, T/, = (i-yi)---(i-r,_i)r, (i ==2,3,...), (i.i) 
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where Fi , F2 , . . . are independent and each Yi is Beta(l — a, + iQ!)-distributed. Let {Vi) 
be the decreasing order statistics of (V^), namely Vi>V2>--- are the ranked values of 
(Vi), and define PD{a,6) to be the law of (V^) on Voo- There is some background to 
the study of these distributions, as explained in [45] and [44]. In particular, a number of 
results concerning PD{a,9) were obtained in [45]. The proofs there use auxiliary random 
variables and related processes (such as stable subordinators and gamma processes) and 
require deep insight into them. As for the original Poisson-Dirichlet distributions, which 
form a one-parameter family {PD(0,6') ■.d>0} and correspond to the gamma processes, 
certain independence property often makes the analysis relatively easier. See [50] and 
[51] for extensive discussions. 

One purpose of this article is to provide another approach to study the two-parameter 
family {PD{a,0) :0 <a < 1,0 > —a}, based mainly on conventional arguments in the 
theory of point processes; see, for example, [7] for general accounts of the theory. This 
means that a random element (Vi) governed by PD(a,0) is studied through the ran- 
dom point process ^ ^Vi , which we call the (two-parameter) Poisson-Dirichlet point 
process with parameters {a, 9), or simply the PD(a,0) process. Note that although the 
above ^ is a 'non-labeled' object, it is sufficient to recover the law of the ranked sequence 
{Vi). For example, we have ^([t,oo)) = precisely when Vi < t. More generally, for each 
n = 2, 3, . . . , quantitative information on Vi , . . . , can be derived from ^ by the principle 
of inclusion-exclusion. Among the many ways of characterizing a point process, we choose 
the one prescribed in terms of correlation functions. For each positive integer n, the nth 
correlation function of a point process is informally defined as the mean density of tuples 
of n distinct points in the process. As far as the one-parameter family {PD(0,6') ■.6>0} 
is concerned, the idea of such an approach is not new. The correlation functions of the 
PD(0,6') process were computed by Wattcrson [54], where these are referred to as the 
multivariate frequency spectra. (See also [39].) 

One advantage of this approach can be described as follows. As Kingman [33] men- 
tioned, the Poisson-Dirichlet distribution is "rather less than user-friendly" . In other 
words, PD(0,6') is not so easy to handle directly. In contrast, the correlation functions 
of the PD(O,0) process obtained in [54] are of a certain simple form. By using them, 
Griffiths [19] obtained several distributional results for PD(0,6'). As a two-parameter 
generalization, we will find ((2.3) below) the correlation functions of the PD{a,9) pro- 
cess. This result will play a key role in revealing the mathematical structure of the 
two-parameter Poisson-Dirichlet distribution via the associated point process. Indeed, 
by exploiting some general tools from the theory of point processes, it will be possible 
to deduce results which extend previously known results for the one-parameter case, for 
example, the joint probability density obtained by Watterson [54] and the moment for- 
mula due to Griffiths [18]. It should be noted that Pitman and Yor [45] essentially found 
the joint probability density in the two-parameter setting. However, we emphasize that 
our expression (5.13) in Theorem 5.4 is quite consistent with Watterson's formula and 
useful for the purpose of further analysis. 

Another aspect of results presented in this article is a development of the Laplace 
transform method discussed in the last section of [45]. Roughly speaking, the main re- 
sult ((3.8) in Theorem 3.2) in this direction states that the (suitably modified) Laplace 
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transform of a 'probability generating functional' of the PD(a,0) process is connected 
with that of a certain Poisson point process on (0,cxj) through some nonlinear function. 
This will provide a powerful tool, especially for the study of certain asymptotic prob- 
lems regarding PD(a, 6) as 9 ^ oo. Such problems have been studied in the case a = 0, 
motivated by the study of population genetics; see [18] for background information. We 
will show that the presence of a fixed parameter a G (0, 1) does not affect the validity of 
results of such type, generalizing a limit theorem of Griffiths [18] and the central limit 
theorem obtained by Joyce, Krone and Kurtz [28]. In the proofs. Theorem 3.2 combined 
with the Laplace method is seen to be very effective in demonstrating such results. 

The aforementioned structure of PD(a, 9) will also be used to explain certain properties 
of the generalized Dirichlet process rj ■=^Vi5xi, where Xi are i.i.d. random variables 
independent of (Vi). The Dirichlet process was originally introduced by Ferguson [17] as a 
prior distribution for nonparametric problems and it corresponds to the case where (Vi) is 
governed by PD(0, 9). In such a context, many attempts have been made to find the exact 
distribution of the (random) mean M ■.= J2^i^i of V- ^ key ingredient for this is the 
so-called Markov-Krein identity, which relates the law of M to the common law v oi Xi. 
Its two-parameter generalization was obtained by Tsilevich [49]; see also [31, 50, 53]. Our 
study of the generalized Dirichlet process will be based on functional analysis of PD(a, 9): 
we regard this process as a 'functional' of PD(a,0) with underlying parameter v. It will 
be seen that Theorem 3.2 explains the mathematical structure behind the Markov-Krein 
identity, a complementary result of which is presented. We would expect that such a 
point of view would help us to discuss further the generalized Dirichlet process. 

The paper is organized as follows. In the next section, we describe the formula for the 
correlation function of the PD{a,9) process. Section 3 deals with the Laplace transform 
method, which will turn out to be another basic tool in the subsequent argument. In 
Section 4, we discuss the Dickman function and its two-parameter version which describes 
the law of Vi with respect to PD{a,9). The main result here ((4.9) in Theorem 4.2) is 
an integral equation satisfied by the generalized Dickman function. Section 5 shows how 
the correlation functions combined with techniques from the theory of point processes 
yield distributional information on PD{a,9). In Section 6, we illustrate the utility of the 
Laplace transform method established in Section 3 by proving some asymptotic results 
on PD{a,9) with 9 large. Section 7 is concerned with the generalized Dirichlet process, 
some of whose properties will reduce to the results obtained in previous sections. 

2. The correlation function of the PD(ck,0) process 

The main subject of this section is the correlation function. We begin with some defi- 
nitions and concepts from the theory of point processes [7]. (See also Section 2 in [27] 
for a comprehensive exposition of this material.) Let ^ be a point process taking values 
in R. For simplicity, suppose that ^ is expressed as ^ = ^ Sxi for some random vari- 
ables Xi,X2, Throughout what follows, we assume that ^ is simple, in the sense that 

Xi ^ Xj{i 7^ j) a.s. For any positive integer n, the nth correlation measure (also called 
the nth factorial moment measure) of ^, if it exists, is defined to be a cr-finite measure 
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such that for non-negative measurable functions / on R" 



E 



f{X,„...,X,J - / fix,,. 



,Xn)^J-n(dXl ■ ■ ■ dXn), (2.1) 



where the sum is taken over ri-tuplcs of distinct indices. In particular, fii is the mean 
measure of ^ and it follows that ^„(dxi • • • d.T„) is necessarily symmetric in xi, . . . ,Xn- 
If /z„ has a density with respect to the n-dimensional Lebesgue measure, the density is 
called the nth correlation function of ^. 

In what follows, a and 6 are such that < a < 1 and 9 > —a, unless otherwise men- 
tioned. Notice that the PD(q;,6') process £, — J2^v ^J2^Vi, with (Vi) and (1^) defined 
as in the Introduction, is simple. Denote by Cn,a^e the product 

fr r(g + i + (^-iM ^ P ' + i)r(0/« + n)a-^ " ~ °' 

In view of the left-hand side, it is clear that 

(m,ne{0,l,2,...}) (2.2) 

with the convention that cq^^.o = 1- We also need the following notation for the n- 
dimensional unit simplex: 

A„ = {(wi, . . .,u„) :ui > 0,. . .,w„ > 0,ui H h Wn < 1}. 

In general, the indicator function of a set (or an event) E is denoted by 1^;. The main 
result of this section generalizes Watterson's formula [54] ((18), page 644) for the corre- 
lation functions of the PD(O,0) process to the PD{a,9) process. 

Theorem 2.1. For each n = 1, 2, . . . , the nth correlation function of the PD(a, 0) process 
is given by 

n / " \ 

qn,a,e{,vi,...,v„) :=c„,Q,e]^w~'^"^^'' 1 ^ 1a„(wi, ■ •■,««)• (2.3) 

i=i V i=i / 

We will give a proof based on a known property of the size-biased permutation for 
PD(a, 9). For the one-parameter family {PD(0, 9):9> 0}, such an idea to derive (2.3) is 
suggested in [39], Corollary 7.4. From the definition of size-biased permutation (see [45] 
or Section 4 of [3]), we can make the following observation. If a sequence [Vi) of random 
variables such that Vi > (i = 1, 2, . . .) and ^ 1^ = 1 is given, then, for each n = 1, 2, . . . , 
the nth correlation measure of ^ exists and 

^„(dwi • • • dvn) = njf'r^ ^1 ^ X]'^'-'^ |ia„(wi, ■ ■■,Vn)tii,{dvi ■ ■ ■ dw„), (2.4) 
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where is the joint law of (V/, . . . , V^), the first n components of the size-biased permu- 
tation (y/) of (Vi). In (2.4) and throughout, we adopt the convention that X]j=i • • • = 0. 

Proof of Theorem 2.1. Let (Vi) have PD{a,9) distribution. A remarkable result ob- 
tained in [37, 41, 43] then tells us that the law of (V^) coincides with the law of (Vi) 
in (1.1). Hence for any non- negative measurable function / on R", the expectation 
E[f{vl...,V^)] is given by 

/ dyi •••d2/„/(z;i,...,z;„)c„.„,en{yz""(l-yO'+*""'}, (2-5) 



where Vi = Vi{yi, . . . ,yn) = Vi 11^=1(1 ~ Vi)^ with the convention that 0^=1 • • ■ = 1- Note 
that the mapping (?/i, . . . , j/„) =: y h- > (ui, . . . , Vn) v from [0, 1]" to A„ has the inverse 



yj = yt(v)=w, l^l-^WjJ (i=:l,...,n). (2.6) 
Therefore, (2.5) becomes 

cn,^fi d«i ••• d«„ 1^ /(v)n{yr(l-y.)''+'""'}- (2.7) 
Observing from (2.6) that 

z=l z=l V j=l / 

we have, by (2.5) and (2.7), 

E[fiVl...,Vi)]= |^^dv^ •••d7;„/(v)g„,„^e(v)n|«.(^l-E".) |- (2-9) 

where qn.a,e is defined by (2.3). With the help of (2.4), this proves Theorem 2.1. □ 

It is known that correlation functions appear in the expansion of the 'probability 
generating function' of a random point process ^ Sxi see Section 5 in [7] for general 
accounts. This functional is defined to be the expectation of an infinite product of the 
form J[[(7(Xi). For the sake of clarity, we provide the following definition of infinite 
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products. Given a sequence {ti} of complex numbers, define 

n 

lim rr(l + i,)e [l,oo], if t,>0 (^ = l,2,...), 

n— >-oo 

n 

lim I f (1 + G C, if the limit exists. 



11(1 + ^0 



In general, it holds that 



OO OO OO 



(2.10) 



Tt—1 ii <■ ■ ■ <?'„ 



n=l h,...,i„{^) 



Also, + admits the corresponding expansion whenever + 1^*1) finite. 

In what follows, Ea^e[ ■ ■] denotes the expectation with respect to PD(q;,6'). 



Corollary 2.2. For any measurable function (p: (0, 1] — > C, 

" OO 

.1=1 

// the above series converges, then 



(2.11) 



EaJ 



OO /. IL If / IL 



1 



(2.12) 



i=i i 



Proof. (2.11) is immediate from (2.10), (2.1) and (2.3) together. Convergence of the 
series in (2.11) allows us to show (2.12) by the dominated convergence theorem. □ 



For instance, if : (0, 1] — > C is a measurable function with support contained in [e, 1] 
for some < e < 1, the series in (2.11) is easily verified to converge. Roughly speaking, 
the assertion of Corollary 2.2 is equivalent to that of Theorem 2.1 (cf. Propositions 2.2 
and 2.3 in [27]). 
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This section is intended to provide a general tool (Lemma 3.1) and to exploit it in the 
study of PD{a,6) processes. It contains a certain inversion formula for Laplace trans- 
forms. Interestingly, in spite of the generality of the formulation, the formula involves 
Dirichlet measures. This seems to show an aspect of the analytic importance of such 
measures. 



Lemma 3.1. Let /3 > and S > — /3. Let R be a function defined at least on a neigh- 
borhood of in which R is expressed as an absolutely convergent series of the form 
R{u) = riu + r2U^ + • • • . 

(i) Suppose that F and G are complex-valued measurable functions on (0, oo) such 
that /g°°ds/-ie-^«^|F(s)| < oo and dzz-(^+i)e-^«^|G(z)| < oo for some Aq > 0. // 
F and G are connected with each other in such a way that 



oo 



AM dss^-^c-^'F{s)^R[\-'^ dz^^G(z) (3.1 



for sufficiently large A > 0, then F{s) coincides with 



(/3 + J>0) (3.2) 



riG(s) + ^r„/ ii— dvi---Avn-i 

,1 = 2 J^r.-! 1 - Vj 

(/3 = = <5) 



(3.3) 



for a.e. s > 0, where the series and integrals converge absolutely for a.e. s > 0. 

(ii) Suppose that G is a complex-valued measurable function on (0,oo) such that 
dz0-('3+i)e-^o^|G(z)| < oo for some Xq > 0. Then the series (3.2) or (3.3) with |G(-)| 
and |r„| in place of G{-) and respectively, converges for a.e. s > and F[s) defined by 
the expression (3.2) or (3.3) satisfies ds s^~^e~'^°''\F[s)\ < oo. Moreover, the relation 
(3.1) holds for sufficiently large A > 0. 

Proof. We consider only the case where /3 + (5 > since the other case can be handled 
quite similarly. To prove assertion (i), we first show that 



dvi ■■■ dvn 
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is finite for a.e. s > 0. Setting, for eacli s > and 7i = 1, 2, . . . , 

A„(s) = {(zi, . . . , z„) : zi > 0, . . . , z„ > 0, zi H h z„ < s}, (3.4) 

we observe, by Fubini's theorem, tliat 



(5+,9n-l 

1 dt;i ■ ■ ■ Avn 

' ^'1 • • • 



rr 1^(^01 / 1 dzi • • • dz„ 



S+Pn-1 



A^^fi zf V U ) 

= / dzi • • • dz„ i^J^ /" dsc"^"! S-^Zj 

Therefore, term- by-term integration yields 

AM ds/-ie-^^Go(s) = ^|r„| A-M dzi^|G(z)| 
Jo \ "'o ^ 

Since the above series converges for sufficiently large A, Go{s) < oo for a.e. s > and thus 
the desired absolute convergence follows. 

Calculations needed for the proof of (i) have almost been done. Indeed, denoting by 
6*1(5) the sum (3.2), one can sec by obvious modification of the above calculations that 
(3.1) with F replaced by Gi holds for sufficiently large A > 0. We thus conclude, by 
virtue of the uniqueness of Laplace transforms, that F{s) = Gi{s) for a.e. s > 0. This 
proves assertion (i). The proof of assertion (ii) is essentially contained in the above and 
is therefore omitted. □ 



At least formally. (3.3) is obtained as the 'degenerate limit' of (3.2), that is, by first 
setting /3 = in (3.2) and then taking the limit as S ^0. Although our inversion formula 
given in Lemma 3.1(i) requires the Laplace transform to have a certain special form, one 
advantage is that it is described in the 'real world', that is, we do not need any complex 
integrals. Considered as prototypes of the inversion formula are integral representations 
of the Dickman function and the Buchstab function (cf. [1], page 22). Both functions 
appeared in asymptotic number theory [4, 12] and are related to PD(0, 1). On the other 
hand, assertion (ii) will be used below to compute the Laplace transform of a probability 
generating function of the PD{a,9) process discussed in [45]. 



1090 



K. Handa 



In such applications, we will employ 
r r(6')(e''"-l), 

Ra.ei^) { r(0 + i){(i - - i}e-\ 



a = Q,0>O, 

0<a< 1,6*7^0, (3.5) 
0<a< 1, 6' = 0, 



where Cq. ~ a/T{l — a). The precise meanings of the power and logarithm in (3.5) as 
complex functions are as follows. Given p e R and t E C\ (— oo, 0], define tJ' = exp(plogi), 
where logt = log|i| +arg<, with argi being chosen in (— tt, tt). As mentioned in Exercise 
1.2.7 of [44], the function Ra^e itself appears in connection with generalized Stirling 
numbers. Note that three expressions in the right-hand side of (3.5) vary continuously 
in (a, 6), even in the limit as — > or a ^ 0. This fact is consistent with the continuity 
of the two-parameter Poisson-Dirichlet family [52], [50]. More importantly, i?a,6» admits 
an expansion of the form 



n t \ V"^ r(6' -t- an)c„ Q,.e „ 
Ra,e[u)=} j —u 



(3.6) 



as long as |m| < l/Ca{~ oo for a = by definition). Combining Corollary 2.2 and Lemma 
3.1, we prove the following formula, a refinement of Corollaries 49 and 50 in the afore- 
mentioned paper [45] by Pitman and Yor. We will use the convention that inf = oo so 
that, for a subset X of R, inf X < oo means X 7^ 0. 

Theorem 3.2. Suppose that g: (0,oo) — > C is a measurable function such that Xa{g) '■= 
inf{A>0: /Q°°dzz^("+i)e-^^|g(z)-l|<oo}<oo. Put 



\:{g) = M{X'>X^{g):^ 



dz- 



j{g{z) [1,00) for all X>X' 



so that, in particular, XQ{g) = Xalg). Then for a.e. s> 0, 



and for X > A* (g), 
X" 



r(0 + i)7o 



1 

TT cxp 



dss" 



Ea.0 



n(i+i.9(.'*^o-ii 



< 00 



.1=1 



dz- 



-(3(^) 



^ a 



; 1-^ dz 



Cq 



a+ 



t(5(^)-1) 



-e/a 



-Xz 



-log 1-^/ dz— ,(.g(z)-l) , 



a = 0, 6'>0, 
0<a<l, 9^0, 
0<a<l, 6 = 0. 



(3.7) 



(3.8) 
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Proof. Consider 



Fo(s):=K 



\{{l + \9{sV,)-l\) 



which admits a series expansion due to (2.11). Noting (3.6), we apply the first half of 
Lemma 3.1(ii) with /? = a, 5 = 9 and R = Ra.e to show that Fq{s) < oo for a.e. s > 0. 
Also, the last half of Lemma 3.1(ii) can be applied to the series expression of F{s) := 
Ea^elYliLi di^^i)] ~ 1 '^^'^^ to (2-12) and we obtain 



A" 



ds i 



9-l„-A 



e-^-^F(s) = R^,e A"" 



dz- 



t(3W-1) 



for sufficiently large A. This extends to all A > A* (g) by analytic continuation, showing 
(3.8). The proof of Theorem 3.2 is thus complete. □ 



This result will be exploited later in a variety of ways by taking 

{l{o.i)('S)i in Section 4, 
exp(— sP), in Section 6.2, 
'0,y(±s), in Section 7, 

where p> Q and ipi^ is the characteristic function of a distribution v. We conclude this 
section with the observation that the probability generating function of the PD(a,0) 
process can also be characterized by an integral equation. 



Proposition 3.3. Suppose that g:{Q,oo) C is as in Theorem 3.2 and set (t){s) = 
g(s) — 1. Then Tlg{s) ■= Ea^g[Y[°Zi g{sVi)] solves the following integral equation for a.e. 
s>0: 



du / dt{u- t)""g(w - t)t'^-^(ng{t) - 1) = 0; 
Jo 



(3.9) 



(i) for a = and 6>0, 

/ng(s) -e [ g{s- t)t^-^ng{t) dt = O; 

(h) for 0<a<l and 9^0, 

dt{s^t)-°'t''Ug{t)- / dtis-ty^t^gis-t) 
Jo 

-a f du f dt{u~t)-'^°'+^^(l)iu~t)t^Ug{t) (3.10) 
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(iii) for 0<a<l and 9 = 0, 

[ dt{s^t)-^Ug{t)- f t-^g{t)dt 
Jo Jo 



(3.11) 

Proof. Given a measurable function / on (0,oo), let / denote the Laplace transform of 
/ and introduce the temporary notation //^(z) = f{z)/z^ for /3 g R. All equalities below 
involving A hold at least for sufficiently large A. 

(i) Set f{s) = s^Ilg{s) so that /i(s) = s^^^I{g{s). We start with a simplified version 
of (3.8) for a = 0, that is, A*/i(A) = T{9) exp(6'0i(A)). Taking the logarithmic derivative 
of both sides with respect to A, we get 

^ ^^^^ -em(^-eg{x) + ' 



A /i(A) ^ A ' A 

and hence /(A) = 9g{\)fi{\). This suffices to prove (3.9). 

(ii) For notational simplicity, let $(s) = ng(s) — 1 and G{s) = Ilg{s). In the case where 
< a < 1 and 0^0, the equality (3.8) reads 

ex' a Co.- 



<i>i_e(A)= l--^0„+i(A) -1 



r{i + e) ' V A° 

Setting f3 = —9{< a < 1), we convert the above identity to 

(r(l - /3)A^ - /3$^+i(A))^/^ = (r(l - a)A" - a^„+i(A))^/". 

By taking the logarithmic derivative, 

$^(A)+r(l-/3)A^-i _ 0«(A)+r(l-a)A"-i 
r(l - /3)A/^ - ^ r(l - a)A" - a0„+i(A) ' 

Here, considering the function l-yiz) :—\/z^ for every 7 < 1, note that the numerator of 
the left-hand side equals $/3(A) + 1;3(A) = G'/3(A) and, similarly, 

^„(A)+r(l-a)A"-i=g„(A). (3.12) 

Consequently, we have 

1„(A)G^(A) - 1^(A)5„(A) - ^^„+i(A)G^(A) + ^5,(A)$0+i(A) = 0. (3.13) 
This proves (3.10). 
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(,,_t)-"tVn<,(t)- / {s-tyt-^Ag{t) 

^" (3.15) 

dt {(s - ty^et"-^ + a{s - ty^°'+^h'^}(f){s - t){ng{t) - 1) = o. 



For < a < 1 and = 0, 

{s-t)-°'dUg{t)-s-''(l>{s)-a f dt(s-0~*"+^V(s-t)ng(t)=0. (3.16) 
Jo 

The proof changes only after muhiplying (3.13) and (3.14) by A and then uses 

/•oo />oo 

l„(A)AG0(A)-i^(A)A5„(A) = l„(A) / c-^M$^(s) - 1^(A) / e-^M<^„(s) 

Jo Jo 

and Ang(A) = 1 + e~^^ dllg{s), respectively. The details are omitted. 

4. The two-parameter generalization of the Dickman 
function 

One of fundamental 'observables' in a point process on R is the 'position of the last 
particle' (cf. [27]), if any. As for the PD{a,6) process, this is nothing but Vi, the first 
component of a PD(q!, (?)-distributed random element (Vi). For special values of a and 
0, the law of Vi was found in various contexts much earlier than Kingman's discovery 
of the Poisson-Dirichlct limit. For example, the Dickman function [12], usually denoted 
p{-), is identified with p{s) ~ ^'o,i(sVi < 1); see, for example. Section III 5.3 of [48] for 
related discussions in asymptotic number theory and Section 1.1 of [1]. Here and in 
what follows, Pa^e is a probability distribution under which (Vi) is PD(a, 0)-distributed. 
Also, a distribution function found independently in [8] and [35] can be identified with 
Pa,o{V^^ < s)- See also [32] ((77), page 14). It is natural to introduce the two-parameter 
version of the Dickman function by 

Pa,e{s) =Pa,eisVi < 1). 



(iii) Differentiating (3.8) with 6^0 and again using (3.12), we get 

l„(A)n,(A) - ig„(A) - |^„+i(A)n,(A) = 0, (3.14) 

which shows (3.11). The proof is thus completed. 

In fact, the generality of the assumptions in Proposition 3.3 makes the resulting equa- 
tions (3.10) and (3.11) rather complicated. Under an additional hypothesis which ensures 
that both g and Ilg are of bounded variation on each finite interval, the following equa- 
tions are derived instead. 

For < a < 1 and 6* ^ 0, 



1094 



K. Handa 



Clearly, Pa,e{s) = 1 for all s < 1. The one-parameter family {po,e : ^ > 0} has been studied 
in connection with both population genetics [18, 54] and the asymptotic theory of the 
symmetric group [22]. For details, we refer the reader to the identity (4.13) in [3] (resp., 
Lemma 4.7 in [1]), where ge{s) (resp., pe{s)) is identical with s^~^ po^e{s) up to some 
multiplicative constant. It also appears in a natural extension [21] of Dickman's result in 
number theory. Another interesting context in which the one-parameter family arises is 
the identification of limit distributions associated with random minimal directed spanning 
trees [40]. 

The aim of this section is to describe consequences for Pa^s which follow from the 
results in Sections 2 and 3. First, a choice of (f> in (2.12) to give an expression for pa g is 
(j){v) = —l[i.oo){sv), with s > being given. With this choice, (2.12) now reads 



p,,,(.)^V ( ^^"'"■"-V n,c,e(g), (4.1) 

ni 

n=0 

where /o,q.6i(') = 1 and where, for n = 1, 2, . . . , 

/„,.,.(.)=/^ n ''';;Vi d«i...d.„. (4.2) 



The reader is cautioned that our notation In,a.e is in conflict with that of [45]. We observe 
that In,a,e{s) = whenever n > s. So, the right-hand side of (4.1) is in fact a finite sum 
taken over n with < n < s. In the case {a,0) — (0, 1), we recover a well-known formula 
for the Dickman function (cf. (1.35) in [1]), while the above identity with 6 — extends 
the formula for PafliV^^ < s) =: Ha{s), obtained in [35] ((3.7), page 730), to all values 
of s > 1. As a two-parameter example, we find in [34] (Theorem 3.3.1, page 164) this 
type of expression for pi/2.i/2: which determines the limit distribution of the maximal 
size of trees associated with a random mapping. 

It is proved in Propositions 19 and 20 of [45] that for all s > 0, 

PcAs) = , \^T^n „^ / ^ -^^ — (4.3) 



r(6l + a)r(l - a) Jo 
and hence the probability density Pa.B{yi G du)/dv is given by 

^-"("^ ■= ne + am-a) ■ ."+1 p..e+.[^^)Mo,i){v). (4.4) 

For later convenience, we remark that underlying (4.3) are 'termwise equalities' shown 
in the next lemma. For every positive integer n and s > 0, let 

V„(s) = {{vi, . . . ,v„) :wi > • • • > w„ > 0,ui H h w„ < s}. 

In particular, V„(l) is simply denoted by V„. 
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Lemma 4.1. For all n ^ 1,2, .. . and s > 0, it holds that 

In,a,e{s)^n dv — In-l.a,e+a[ • (4.5) 

Jmin{s-\1} + \ V J 

Proof. Obviously, we may assume that s > 1. Also, (4.5) with n = 1 is clear from the 
definition (4.2). For n > 2, by symmetry of the integrand in (4.2), In,a.e{s) equals 

d^n , f dz^i— du^ T \ iAa\ 

iv„-i(l-«„) ■■■Vn-1 \ J^l ) 

The change of variables = (1 — u„)ui converts the inner integral to 



n-1 \0+an-l 
U-i 



1 /I \6l+a-lr '1 



(ri-l)! 

Substituting this into (4.6), we get (4.5). □ 

It should be noted that (4.3) does not give a closed equation for pafi unless a = 0. To 
derive such equations for the general case, we shall apply some results from the previous 
section. Setting g = l(o,i) in (3.8) and observing that A* (5) = 0, we see from Theorem 
3.2 that for aU A > 0, 

\' dss'-^c-^'{p^,e{s)-l) = R,,^g(-X-" dz^y (4.7) 

For 9 > 0, this relation is rewritten in a slightly simpler form: 



r(^) 



exp(-(?/rdz^), 



Q! = 0, 



\e poo 

ds/-ie-^V„,e(s)=<( ^ ^ ;_e/a (4.8) 



l + ^/rd;.f^) , 0<a<l. 



Integral equations satisfied by p^.g will be derived as the ones equivalent to (4.7), in 
other words, as consequences of Proposition 3.3, its variant (3.15) and (3.16). Preliminary 
observations concerning this are that ng(s) = Pa.d{s) for g = l(o.i) and that (3.9), (3.15) 
and (3.16) hold for all s > 1 because of continuity. As will be discussed, the equation 
below is known in the cases a = and 9 = 0. For this reason, the proof of the next 
theorem will be concerned with the remaining case only. 
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Theorem 4.2. pa_e solves the equation 

(s-l)V,e(s-l)+ / (s-0""i'dp„,e(t)=0 (s>l). (4.9) 

Js-l 

Proof. The case where < a < 1 and 9^0. (3.15) now reads, for all s > 1, as 

r{s-t)-^tUp^,e{t) + {s-lf 
Jo 

+ r \t{{s-tyet'>-^ + a{s-t)-^''+^h<'}{pa,,e{t)-l)=0 
Jo 

because —dg(t) on (0,oo) equals Si{dt). Integration by parts yields (4.9). □ 

Remarks, (i) (4.9) with a = or its variant can be found in [19, 21, 22, 54]; see also 
(4.25) in [1]. 

(ii) One can easily see that (4.9) with 6* = is derived from a functional equation for 
Ha{t) = 1 — Pa,o{t) obtained in [35] ((3.5), page 730). 

In the next section, we calculate not only the marginal distributions Pa.eiVm G •) for 
m = 2, 3, . . . , but also multidimensional distributions of PD{a,0), by developing point 
process calculus based on Theorem 2.1. 



5. Distributional results for PD(q;,0) 

In this section, we apply the theory of point processes to deduce from (2.3) some distribu- 
tional information on PD{a,d). An essential idea underlying the subsequent argument is 
the principle of inclusion-exclusion, which was already being used in [19]. By the same rea- 
soning as in the proof of Theorem 2 there, namely, by a version of Frechet's formula (see, 
for example. Section IV. 5 in [15]), the following relationship between onc-dimcnsional 
distributions and correlation functions holds with great generality. 

Lemma 5.1. Let — oo < a < 5 < oo and let ^ be a simple point process on (a, 6) with 
correlation measures /ii,/^2,---- Suppose that ^((a,&)) = oo a.s. and that ^([c,6)) < oo 
a.s. for each c£ {a,b). Given a positive integer m, let Zm be the mth largest point in ^. 
Then, for every z € (a, 6), 

P{Z„i >z) = — - V — — / p„i+k{<iyi ■ ■ ■ dym+k)- (5.1) 

{m - 1)1 ^ (m + k)kl 7[2.f,)™+fc 

//, in addition, ^ has correlation functions (71,(72,..., then 

P(z,„edz) _ i(„^b)(^)>^(-i) 



dz (m - 

^ ' k=0 



'z) °° {-lY f 

7TT ,1 / dyi ■■■ dy.m+k-lqm+k{z,yl,...,ym+k-l)■ 

^)■ t^o Jlz,b)'^+''-' 
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Example 5.1. Consider a Poisson point process Zi > Z2 > ■ ■ ■ on (a,&) with mean 
measure A(d2:) = h{z)dz such that A((a,6)) ~ 00 and A([c,6)) < 00 for aU cS {a,h). 
Lemma 5.1 then gives 

P{z^ e dz) h{z) ' ^ "-^ 



dz (m - 



hiy)dyj exp(^-^ /i(y)dy)l(,,6)(z) (5.2) 



because the nth correlation function is h{yi) ■ ■ ■ h{yn); see, for example, Example 2.5 in 
[27]. 

Setting Pm,a.e{s) ~ Pa.e{sVm < 1) {m = 1,2,...), we obtain the fohowing proposition 
containing a two-paramctcr generalization of the aforementioned result in [19]. 

Proposition 5.2. Let m be a positive integer. Then, for all s > 0, 

1 

Pm,a,e(s) = 1 — 7 TTT 7 , , , , Cm+k,a.9 Im+k,a.d (s) (5.3) 

(m—iy. ^-^ {m + k)k\ 



and for all A > 0, 

X' dss^-ie-^^(p„,„,e(s)-l) = i?„,„,e(^-A-"y^ dz^], (5.4) 

where 



D ^ \ \.J I'm— 1,0,9 / m — 1 D N J 

(m-1)! Jo 



r(. + i)(-.)™- r,™-i,..,,, 



(5.5) 



("^-1)! Jo 



■ dx, < a < 1. 



(m - l)!r(e'/a + l)r(l - a)™ 7o (1 - Caa;)''/°+^ 
^/so, for V e (0, 1), 

du r(6l + a)r(l -a) ^ 

(5.6) 

1 (-1)'' , fl-v' 

X7 TTT / , Ti C„i+fc-l,a,e+QJm+fc-l,Q,e+Q 

m — 1 ! fc! ' . I 

^ ' 0<fe<l/t>-m 

Remarks, (i) Since pi^a.e = Pa.e, it is worth noting that the right-hand sides of (5.3), 
(5.5) and (5.6) with m = 1 equal the right-hand sides of (3.5), (4.1) and (4.4), respectively. 

(ii) In [34], Theorem 3 on page 47 gives the expression (5.3) with [a, 9) = (0, 1) to the 
limit distribution of the mth maximal cycle length in a random permutation. 



1098 



K. Handa 



(iii) At the end of this section, another expression of the density (5.6) wiU be given in 
terms of the generahzcd Dickman function. 

Proof of Proposition 5.2. First, (5.1) and Theorem 2.1 together imply that 

1 °° (-1)*"' 

Pa,eisVm > 1) = 7 TTT > , 7 , , n , , Cm+k,a.e Im+k,a.9 (s) ■ 

(m — 1)! ^— ^ (m + k)k! 

By removing the terms which actually vanish, (5.3) follows. Next, showing (5.4) reduces 
to verifying that for u sufficiently close to 0, 

R . , (-l)"'-^ f> r(g + a(m + fc)) , 

For, if (5.7) is true, then (5.3) and Lemma 3.1(ii) with R ~ Rm,a,e and G = — l[i,oo) show 
that (5.4) is valid, at least for A sufficiently large, and the expression (5.5) allows us to 
extend (5.4) to all A > by analytic continuation. Turning to (5.7), one can easily verify 
it by using (term-by-term differentiation of) (3.6) or, alternatively, by substituting the 
expansion of e^^ or (1 — Cq2:)~(^/"+™^ into (5.5). This proves (5.4). 

Lastly, since (5.3) is at hand, the proof of (5.6) is similar to that of (4.4) based on 
(4.5). The details are left to the reader. □ 

To obtain multidimensional results for PD{a,6), one needs more developed arguments 
still based essentially on inclusion-exclusion. In the theory of point processes, one such 
calculus is formulated as a connection between correlation measures and Jonassy mea- 
sures; see Section 5 of [7]. Its significance is that the local probabilistic structure of points 
in the process is revealed in terms of correlation measures. 

Lemma 5.3. Let — oo < a < 5 < oo and let he a simple point process on {a,b) with 
correlation functions qi,q2,.... Suppose that {{a, b)) = oo a.s. and that for cached {oL,b), 
the nth factorial moments Mn{c) o/^([c, 6)) satisfy 

f^^!^^{l+eY <oo for some e = e{c) > 0. (5.8) 

n=l 

Let Z\> Zi> ■ ■ ■ he the decreasing sequence of points in ^. Then, for each m— 1,2,..., 
the joint probability density of {Zi, . . . , Z^) is given by 

/„i(z)=V— — / dyi ••• d?;fc(7m+fc(z,2/i,...,2/fc), (5.9) 



fe=0 



where z = (zi, . . . , z„i) is such that b > zi > ■ ■ ■ > z„i > a. 
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Proof. Let 6 = : zq > zi > • • • > Zm > a be given arbitrarily. Consider the localized process 
^("^(dx) :=^(dxn [z„i,6)), the nth correlation function of which is given by gi™'(x„) := 
g„(x„)l[^^^,,)(a;i) • • • where x„ = (xi, . . . ,a;„). By applying Theorem 5.4.II in 

[7], then, we have 

P{Zi > Zi > Z2 > Z2 > ■ ■ ■ > Zm-2 > Zm-l > Zm-1, Zm > Zm) 

= P(^(™)([z,„,Zo]) =7n- l,C("n(^l,^0]) = • •• = C^™H(Z„,_1,Z,„_2]) = 1) 

00 / ^s^i^ . m— 1 

= ^^T~ / n '^{z,.z,-l]i^^hln-l+ki^rn-l+k)dxi ■■■ dx,n-l+k 

dx,n~i ■ ■ ■ dx2 / dxi 

Zm-1 ^ J zi 



71— dXmdXm+1 ■ ■ ■ dXm-l+kqm-l + k{^m-l+k) (5.10) 



oo 

X 

k=0 

Zm-2 



d.T,„_i • • • / dx2 I dxi 



^ nb pb pb 

X ^(-1)'' / dxm j dxm+i---j dxm„i+feg„_i+/c(x,„_i+fc), (5.11) 



where the last equality is due to symmetry in Xm, ■ ■ ■ ,Xm-i+k and the terms for 
fc = in (5.10) and (5.11) arc imderstood as qm-ii'^m-i) with the convention that 
go = 1. This shows (5.9) by the symmetry of the integrand in (5.11) in A: — 1 variables 

Example 5. 2. If ^ is a Poisson point process with mean measure density h as in Example 
5.1, the joint density of the first m largest points in ^ is given by 

E^i^ / dyi---dykl[h{z,)l[h{y,)=Y[h{z,)cxp(- h{z)dz), (5.12) 

fe=0 ''[zm,b)- \ Jzm / 

provided that 6 > zi > • • • > z,„ > a. 

By applying (5.9), we deduce the next result for joint densities, which is essentially 
contained in [45] , although an explicit formula is not given. (Indeed, the authors obtained 
in their Corollary 41 a corresponding result for the variables ?/i,j/2,-- - in (2.6) instead 
of i'i,W2,.... Alternatively, the formula below can be retrieved from their Proposition 
47 after some additional calculations. See also Lemma 3.1 in [16].) This simultaneously 
generalizes the formula due to Watterson [54] for a = and the one-dimensional result 
(4.4). Also, our expression (5.13) below will be quite useful in Section 6.1. 



1100 



K. Handa 



Theorem 5.4. Let (Vi) be governed by PD(a, 9). Then, for each m = 1, 2, . . . , the joint 
probability density of (Vi, . . . , Vm) at v = (vi, . . . , w,„) € Vm is 

/m,a,e(v) := C™,a,e 1 - Pa.e+Qm ( '^"^ 1 ■ (5.13) 

i=l \ j = l / \ m / 

Proof. First note that the PD(a, 9) proecss £,='}2i5vi is a point proeess on (0, 1) which 
satisfies aU of the assumptions in Lemma 5.3 since, for each < c < 1, ^([c, 1)) < 1/c 
a.s. by Vi<l/i (i = 1, 2, . . .). According to (5.9), the density to be computed is 



(—1) [ 

— -rr- I dui ••• dwfcg„+fc,a,e(v,Mi,...,Ufe). 

Here, by (2.3) and (2.2), 

dui ••• dMfcg,„+fc,a.e(v,'Ui,...,Mfc) 

(5.14) 



i=l 

k ^ / m X e+a(m+fc)-l 



with notation (3.4). Because the last integral in (5.14) is equal to 

6+am— 1 



the desired expression (5.13) is derived from (5.14), (5.15) and (4.1). □ 

Since fm,a,e is too complicated to compute quantities concerning the joint distribution 
directly from this, it is worth providing a moment formula. For the one-parameter case 
where a = and 6* > 0, Griffiths [18] showed that Eo^elVi"^ ■ ■ ■ Kn""] equals 



f dzf[{z^'-\-'^)cxp(-9 f dz—), (5.16) 



r(6' + fli H h a„i) y,^>...>^^>o 

provided that 9 + ai + ■ ■ ■ + Om > 0, where dz = dzi • • • dzm- The proof reduces to calculus 
of a gamma process by virtue of the well-known independence property of it; see [50] and 
[51] for extensive discussions and related topics. In this regard, we remark only that the 
integrand in (5.16) contains the density function, computed from (5.12), of the first m 
largest points in a gamma process with parameter 9, that is, a Poisson point process on 
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(0,00) with mean measure 9dz/{ze^). Although the case where < a < 1 and 9 > —a 
can be handled by means of representing a PD(a, 0)-distributed random element by an 
Q-stable subordinator ([45, 50] combined with (5.12)), we prefer to exploit the previous 
results in order to make the proof self-contained. 

Corollary 5.5. Let < a < 1 and 6 > —a. // oi H + a™ > —9, then 

r(^ + i)r(0/a + m)a™-i ^^^^^^ 



r(6' + ai + • • • + a,n)T{9/a + 1)T{1 - a)™ 

-(9/a+m) 



X 



r '"■ / roc ~z 

/ dzn(^r"^"^^'o-) 1+a / dziL^ 



Proof. By (5.13), 



= / dui • • • dt;™ ]^ Pa.e+ami j- 

■^V™ i=i V 3 = 1 J \ J 

Multiplying both sides by T{9 + A) = y^+^~^e~y dy with A = ai-\ h a^, we intro- 
duce an additional integration with respect to a new variable y on the right-hand side. 
For the resulting (m -I- l)-dimensional integral, perform the change of variables 

s :=Wrn^ ^1 - XJ^J^ ' zr-=v,y (i = l,...,m), 

or, equivalently, y = (zi H h z„i) + z,„s, Vi = Zi/y (i = 1, . . . , m), to get 

E^m[Vi''' ■ ■ ■ Vm.^"^]r{9 + A)/Cm.^^^e 

= dzn^r"'"''"^^"""' / dsy™e-«|J|/+— Va,.-Ka™(s) 

J2i>->2,„>0 j^j^ Jo 

with Jacobian J — d{vi, . . . , Vni,y)/d{zi, . . . , Zm, s) = y^™Zm- Thanks to (4.8), with 9 + 
am > in place of 6*, we arrive at (5.17). □ 

Remarks, (i) It can be seen that the formula for Eafi\V^^-]{p > —9) in Proposition 17 of 
[45] is recovered by putting ai ~ ■ ■ ■ = a„j_i = and a„j = p in (5.17). 

(ii) Let < w < 1 and m € {2, 3, . . .}. By integrating /m,Q,e(wi, . . . , Vm-i,v) with respect 
to dvi ■ ■ ■ dvm-i over {1 > ui > • • • > Vm-i > v}, we can deduce a representation (an 
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alternative to (5.6)) of the density Pa,e{Vm € dv)/dv of the form 



where m„i = 1 — (ui + • • • + u„i_i). This formula, a natm'al extension of (4.4), is verified 
by calculations similar to (5.14) and (5.15). The details are left to the reader. 



6. Asymptotics of PD(ck,0) for large 6 

In this section, we study certain asymptotic behaviors of PD(q!,0) as — > oo, gener- 
alizing results of Griffiths [18] and of Joyce, Krone and Kurtz [28], who all worked on 
PD(0, 6')'s with motivation coming from the study of population genetics. Although there 
is a context [24] in which such an extension could be applicable, this section is mainly 
intended to demonstrate that the results we have thus far obtained provide efficient meth- 
ods for the study of two-parameter Poisson-Dirichlct distributions. It turns out that the 
presence of the parameter a does not affect the validity of the assertions analogous to 
those for PD(0, 0)'s, except for some minor changes of sub- leading terms or multiplicative 
constants in rescaling. 



6.1. Convergence to a Gumbel point process 

Given < a < 1 and 6 >1, put 

Pc.,e = \oge - (a + 1) logloge - logr(l - a). (6.1) 

Under the assumption that (l^^^^)^i is distributed according to PD(O,0), it was shown 
in [18] that as 6* — >■ oo, [QVl ' — f3o,e)^i converges in law to a Poisson point process with 
mean measure density exp(— z), — oo < z < oo, which may be called a Gumbel point pro- 
cess since the largest point in the process obeys the Gumbel distribution cxp(— cxp{—z)), 
as seen from (5.2); see also (6.3) below. A two-parameter generalization of this result is 
the following. 

Theorem 6.1. Fix an a e [0, 1) arbitrarily. For each >l, let (V^^"'^^)^! have PD(a, 9) 
distribution and define 

^(<.fi)^Qy(<.,9) (i^l,2,...). (6.2) 

Then, as 9 ^ oo, {z[°''^\ Zj"'^', . . .) converges in joint distribution to the decreasing order 
statistics {Z^, , . . .) of a Poisson point process on the whole real axis with mean measure 
exp(— z) dz. 
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The proof of this theorem reduces to showing pointwise convergence of joint densities, 
thanks to SchcfFc's theorem (cf. [2]). It follows from (5.12) that 



P{Zl edzi,...,Z^e dzm) = exp {~y^Ji 



dzi ■ ■ ■ dzn 



(6.3) 



i=l 



while the density of {z[°''^\ . . . , Zm'^^) derived from (5.13) involves the generalized Dick- 
man function as a key factor. It is therefore essential to study the asymptotic behavior 
of pa,e with a properly rescaled variable for large 0. This is the content of the following 
lemma, in which is derived the Gumbel distribution, the law of Z^. 



Lemma 6.2. Fix < a < 1. Let 6*0 > 1 and b: [6'o,oo) R 6e bounded. If V{ 



and 



zj"'^' are as in Theorem 6.1, then for each x € R 



x + f3o 



-HO) 



(6.4) 



as —> oo. In particular, zj"'^^ — > Zl weakly and Pafi{s) — > 1 for all s G R. 

Proof. First, the equality in (6.4) is clear from (6.2). Rather than using (4.1), which is 
not very informative for the estimation of the value of Pa,e^ we employ (4.8), that is. 



Jo 



ds i 



-l-\s 



Pa,e{s) = 



r(6')exp -e dz , 



T{e) 1 + a 



ya+l 



-e/a 



a = 0, 



< a < 1, 



(6.5) 



where A > is arbitrary. We will give only a proof for the case where < a < 1 because 
the case a = can be handled in much the same way. By Fubini's theorem, the left-hand 
side of (6.5) is written as 



dss'^-^c-^' 



- En 



dtt'-^c-' 



EnAhAVi-') 



where fx,eis) = J^'^ dtt^ ^e *. For all A > and s > 0, we get 

fxAs) .fxAs) V Jx 

as a Chebyshev-type bound. Set e{9) = 6^^^^, so that 

£(6i)iog6'->o, e{e)Ve 



(6.6) 



(6.7) 
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Given arbitrary a; G R, choose 

s = s{x,e):= ^-—+bie) and A = A(a;, 6*) (1 + e(6'))(x + /3„ e), 

both of which are positive for sufficiently large 9. Then £{x,9), defined implicitly by 
X{x,9)s{x,9) = {l + e{x,9))9, has the same properties as (6.7). The latter property, com- 
bined with the standard argument in the proof of Stirling's formula (that is, the Laplace 
method [9]), shows that 



as 9 ^ oo. Consequently, by (6.6), with the above choice of s and A, we have 



limsuppQ,e — — h &(6') < limsup 1 + Cq / dz ——r ] . (6.8) 

For two functions ai{9) and a2{9), which may have parameters a, x, etcetera, the notation 
01(6*) ~ 02(6*) wiU mean that ai{9)/a2i9) ^- 1 as 6^ ^ 00. By (6.7) and (6.1), 

g-z e-A(x,e) g-(:r+/3„,«) r(l-a) _^ 
M^,e) ^I^^A(x,0)"+i"(a: + ^„,e)"+i ~9 ' 

Combining this with (6.8) yields 

limsupp^.e ( — ^ + bi9)] < exp(-c-"). (6.9) 

The converse estimate can be shown in an analogous way. Indeed, considering this time 
1 — Pa,e{s) = Pa,e{Vi > s"""^), we again have, by Fubini's theorem, 

/>oo 

A"/ dss'-^e-^'{l-p^,g{s))=E^,g[<^x.e{Vi)l (6.10) 
Jo 

where '^\,e{s) — di i^~^e~*. Therefore, a Chebyshev-type bound we get is 

._.„,.,>,^^._^|._(..c./;.^)-- 

where the last equality follows from (6.10) and (6.5) together. By setting s = s{x,9) and 
A = (1 — £{9)){x + Pa,9), one can easily modify the previous argument to obtain 



limsup-^ 1 - pa^e h b{9) ]><l- cxp(-c"'=). 
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This, together with (6.9), proves the convergence in (6.4). The other assertions can be 
seen simply by taking b{-) = and noting that Pa,ei') is non-increasing. □ 

Proof of Theorem 6.1. Fix a positive integer m and zi > ■ ■ ■ > z,n arbitrarily. By 
(5.13), the joint probability density of {z[°''^\ . . . , Zm'^^) at (zi, . . . , Zm) is 



Jrr 



/ m,a.l 



^1 ^" l^a^O Zyji 



\ e+am-1 



- C,n,.fi Wiz^ + /3o..)^'"+^^ 1 - E ^^^^ (6-11) 

It is easy to see that the factors in (6.11) behave as 

m 

C™,o,en(^»+/3a.e)"'"''"^^^*'""''"r(l-«)^™(log0)-(l+")™ 
1=1 

and 



(^1 - E j ^ exp ^- E j 0-"(log0)(i+")"r(l - ar, 

respectively. In view of (6.3), it remains to show that the factor involving Pa,e+am{' ■ ■) 
converges to exp(— exp(— 2;^)). But this follows from (6.4) with &(•) defined appropriately. 
The proof of Theorem 6.1 is thus complete. □ 

The weak convergence of Z^'^' to implies the following. 

Corollary 6.3. For each m = 1, 2, . . . and .t e R, as 9 ^ 00, 

P_.(^)-oxp(-e-^)g£^. (6.12) 

Proof. As a consequence of Theorem 6.1, the left-hand side of (6.12), that is, Pizln'"^ < 
x), tends to P{Z*^ < x), which is calculated by (5.2): 

P{Z:, < x) = r exp(-my - e-^) = H t—ie"* di. 

J-oo[Tn-l)\ (m-l)!7„-x 

This changes into the right-hand side of (6.12) after repeated integration by parts. □ 
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Theorem 6.1 immediately implies that (V/"' , ,...)—>■ (0, 0, .. .) as 6* — > oo. Feng 

[16] has recently obtained associated large deviation estimates with rate function {vi) i— >■ 
— log(l — '^Vi). The proof is based on logarithmic asymptotics of the finite-dimensional 
densities. We note that, in view of (5.13), one can grasp the validity of such a result 
because the factor involving Pa,e+am stays away from by (6.4). 

6.2. Asymptotic normality of the population moments 

As the final topic on PD{a,9) for large 9, we discuss asymptotic normality of the popu- 
lation moment, which is defined for each p > by the sum 



where (V\ ,V"2 , . ■ •) continues to be PD(a, 0)-distributed. We can verify that 
oo, a.s. for all p> a. Indeed, setting gp{s) — exp(— s^"), one observes that Xa{gp) 
in Theorem 3.2 is finite for p > a and the claim follows from (3.7) with g = gp. Also, 
(2.12) with (j){v) = gp{sv) — 1 makes it possible to derive the Laplace transform of the 



law of Hp' because Y^Lidpi^^i' ) = exp(— s^'i/p ' ). As far as the PD(O,0) case 



is concerned, the statistics H2 ' , called the population homozygosity, is known to play 
a special role in some population genetics contexts (see [19, 20, 29]), where V^^'^^ are 
regarded as the ranked frequencies of alleles in a population. More generally, Hp^'^^ is re- 
ferred to as the pth population moment in [28], where the asymptotic normality of Hp^'^^ 
for large 9 has been established. On the other hand, covariances of Hp°''^^ 's are calculated 
in [10] ((12), page 188), motivated by certain physical problems. In this subsection, we 
will provide an extension of these results to the two-parameter setting. 

In the case a = 0, it was shown in [28] that for each p = 2,3,..., as 6*— >oo, the limit 
distribution of Wp^^ := \/9{9p^^ Hp^'^^ /T{p) — 1) is the normal distribution with mean 
and variance r(2p)/r(p)^ -^P^- In fact, in [28], much stronger results were obtained, 
showing, for example, asymptotic normality of random vectors {Wm )m=2- To seek an 
appropriate rescaling for the two-parameter case, we give some auxiliary results by means 
of the first and second correlation functions, namely, qi^a,e and q2.a,e- 

Lemma 6.4. Let p> a and p' > a. Then 




r{9 + i)r{p-a) 

T{9+p)r{l-a) 




(6.13) 



and the covariance Cov{Hp' \H^' ) of Hj 



(a,e) 



and H[ 



r(a,e) 

p' 



is given by 



p 




(6.14) 
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Moreover, as 6 ^ oo, 



and 



^ P ^ r(l-a) 



(6.15) 



Co.iHr\K?^^^) ^ + (6.16) 

Proof. (6.13) and (6.14) arc immediate by calculating /J^ f^<7i,Q.e(v) du and 



respectively. By virtue of the well-known formula r{d + l)/r{9 + p) ^ O^^p implied 
by Stirling's formula, (6.15) follows from (6.13). Showing (6.16), however, requires a 
more accurate version of Stirling's formula: with the standard o-notation, T(0 + 1) = 
6^\/2n9exp{—d + C6^^ + o{9^^)) as 6* — > oo for some universal constant C; see, for ex- 
ample, [9]. This makes it possible to show that 

^J|±i| = e'-p cM-pip mer' + o{9-')). (6.17) 

The proof of (6.16) is a matter of straightforward calculation, combining (6.17) with 
(6.14). The details are omitted. □ 



Lemma 6.4 suggests a rescaling of the form 

(6.18) 

We can now state a desired result for < a < 1 . 



P \T{p-a) P 



Theorem 6.5. Fix an a € (0,1) arbitrarily. Suppose that p> a and p^l. For each 
> 0, let Wp"'^"^ be defined as above. Then, as — > oo, Wp°''^^ converges in law to a 
normal random variable with mean and variance 

o r(l-a)r(2p-a) n „ 

1 (p — ay 

Proof. Fixing an arbitrary x G R, we will prove that as oo, 
V'^"'^)(x) :=^;[exp(v^xWp("'«))] ^exp(^-^ 



P ^2 



cc^ . (6.20) 
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The proof is based on (3.8), a consequence of Theorem 3.2 which can be written as 



r(^) Jo 



.1=1 



dz^(l-g(z)) , (6.21) 



provided that 1 > A* (g). Let g be of the form g{z) — exp(cizP + diz) and observe that 



.i=l 



:S[exp(cisfi/("''')+dis)] 



(6.22) 



Note that (6.21) is vahd as long as the integral on the right-hand side is small enough in 
modulus. With the choices 



^ X r{l- a) ^ X 

Ci = v — 1— — -, r and a\ = —\J — \—;=, 

r{p-a) ^ V^' 



(6.23) 



this requirement for the validity of (6.21) is fulfilled for 9 large enough since as 6* — oo, 

Td-a) „ 



dz- 



y/eJo 0"+iVr(p-a) 



r{p - a) ' 



zP-z] +0(61-1) 



261 



r(l-«)(a„,/ + (p- 1)2) + 0(0-1), 



where the last equality is achieved by direct calculations. This also implies that the 
right-hand side of (6.21) tends to exp{—aa,p^x^/2— {p — \Yx'^/2). On the other hand, 
the left-hand side of (6.21) with (6.23) becomes, by (6.22), after change of variable s = Ot, 



du«-ic-'^*£;[exp(ci6if i/^"'^)tf + diet)]{T{e)e-'>)-^ 

dtt-ie-'"'(*V^"'«)(a;tJ')exp(v^a;V^(tP-<))(r(6l)6'-''e^)-\ 
where h{t) = t — \ogt — 1. Thus, by Stirling's formula, (6.20) is equivalent to 

dtt-^e~^'^W^l-^'\xtP) eMV^-VO{t' - t)) (^.^-iP-if.V2 ^ 

the proof of which will be divided into proofs of 

/>oo 

d^^-le-'"'(*'(V'l,"'''^a;^'')-V'^"'^'(a;))exp(^/=Ta;V^(^P-^))^0 (6.24) 
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d^^-lc-«''(*)cxp(^/^x•V0(^^-^))-c-(p-l)'-'/2_^o. (6.25) 

V 27t7o 

Note that the integral in (6.24) ean be bounded in modulus by 

) />oo 

dtt-ie-'"'(*'|V^"^'''(xtP)- ^^,"'"^2^)1 < \x\E[\Wj,°'''^^\] / d«-ie-^^(*^|tP-l| 

Jo 

and that sup^^j^ < oo sinec, by Lemma 6.4 and (6.17), 

< Var(W^("'«)) + {E[Wl"^o^]f ^ a^^. 

Therefore, our task is reduced to establishing that Vd dtt~^c~^'^''*^tP — 1| and 
(6.25). Both are purely analytic problems which can be solved without difficulty by the 
Laplace method (or by the change of variable t= 1 + u/\/9). 

It remains to show that ctq.p^ > 0. A stronger result which can be shown is 



<^a,p^ = (1 ^ Q;) Var 



r(i-«) 
r(p - a) 



yP- 



1-a 



(6.26) 



where Ya has the density "e ^/r(2 — a), < z < oo. The proof of (6.26) is fairly direct 
and therefore omitted. The proof of Theorem 6.5 is thus complete. □ 

Remark. The above proof also works in the case a = if (6.21) is replaced by the 
corresponding equality for a = 0. This gives a proof which is completely different from 
that in [28] and which is not based on the independence property underlying PD(0,6'). 

One may consider an extension of Theorem 1 in [28] in the two-parameter setting, 
that is, convergence of {Wp"'^'' :p > a} to a centered Gaussian system as — s- oo. The 
convergence of every finite-dimensional law is shown by a suitable modification of the 
proof of Theorem 6.5. It means, for instance, that g in (6.21) is taken to be 



g{s) = exp 



T " 



k=l 



r(i-a) ^ 

T{ph - a)' 



where ife S R, pi > • • • > p„ > a and n = 1, 2, . . . are arbitrary. The limit process is 
described by its covariance V{1 — a)T{p+p' — a) / {T {p — a)T {p' — a))+a — pp' =: C{p,p') 
and the associated quadratic form J2C^{Pk,Pi)xkXi has two different expressions: 



dz- 



r(i-a 

= (l-a)Var(^a:fc 



fT{l-a) 



r(i-a) 



r{pk - a) 



\T{pk - a)' 



yPk 



l-a 



'^XkiPk - 1) f > 
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where is as in (6.26). We do not produce the whole proof, a routine matter which is 
left to the reader. 

7. The generalized Dirichlet process 

This section discusses the two-parameter generalization of Dirichlet processes from the 
point of view of our previous results. There are a number of motivations to study the 
original Dirichlet process, for example, as a prior distribution in Bayesian nonparametric 
statistics [17] and as a stationary state of a certain diffusion process arising in population 
genetics. The well-known relationship between the Dirichlet process and the Poisson- 
Dirichlet distribution is described as follows. Let {Xi}°^^ be a sequence of i.i.d. random 
variables and {Vi)°^i be governed by PD(0,6'). Assuming independence of {Xi}°^^ and 
(^i)i^ii ^6 know that a random distribution 77 :=^V^;(5xi defines a Dirichlet process 
with underlying parameter measure Ov, where ly is the common law of the X^'s. In this 
sense, PD(0,6') is the simplicial part of a Dirichlet process. By replacing PD(0,6') by 
PD(a,0), a two-parameter generalization of a Dirichlet process was introduced in [49]. 

In the case of Dirichlet processes, there are many articles, including [5, 11, 23, 36, 
46], in which exact forms (the distribution function, for example) of the law of the 
(random) mean "^XiVi of 77 are obtained in terms of i^. A key tool is an integral identity 
due to Cifarelli and Regazzini [5] which connects these two laws. As explained in [31, 
50; 53], it is called the Markov-Krein identity because in the case 9 = 1, an integral 
transform is involved, analogous to the one studied by Markov and Krcin in the context 
of moment problems; see [30] for background and various applications. Its extension 
was proven by Tsilcvich [49] for the generalized Dirichlet process; see also [31] for a 
further extension and [50, 53] for a simple proof under some restriction on the support 
of V. This identity gives a one-to-one correspondence between v and the law of '^XiVi 
with (Vi) having PD(a,0) distribution. However, such a correspondence is implicit and 
seems very subtle in general. One therefore needs some procedure of inversion in order 
to obtain explicit information. Recently, James, Lijoi, and Priinster [25] obtained some 
distributional results for the generalized Dirichlet process by means of a Perron-Stieltjes- 
typc inversion formula. Since PD(0, 9) and PD(a, 0) are related to gamma processes and 
stable subordinators, respectively [45], the corresponding problem is naturally considered 
for a more general class of random distributions derived from subordinators. as posed in 
[53]. For results of this kind, we refer the reader to works by Regazzini, Lijoi and Priinster 
[47] and by Nieto-Barajas, Priinster and Walker [38], both of which make essential use 
of the Gurland inversion formula. 

Our focus now will be on what is implied by Theorem 3.2 in the aforementioned context. 
It will be pointed out that the basic identity (3.8) involving the probability generating 
functional exhibits a mathematical structure underlying these kinds of identities. This 
should be compared with the proofs in [11, 31, 49], where the Markov-Krein identity is 
viewed as a relation between moment sequences through the Ewens (-Pitman) sampling 
formula [13, 42]. Let us introduce some notation used to describe the domain of a map 
defined via the Markov-Krein identity. Denoting by V the totality of Borel probability 
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measures on R, let 



Vo^^veV:j i/(da;)log(l + |x|)<oo| 



and set Va = {v : /j^ i^(da;)|x|" < 00} for < a < 1. The Markov-Krein identity al- 
ready mentioned is the following. Taking a.v& Va, let {Xi\°^^ be i.i.d. random variables 
with each law of Xi being v and let (V^i)iSi be PD(a, 0)-distributed. If {Xi]°^^ and 
{Vi)°^i are mutually independent, then ^ < 00 a.s. (as was shown in [14] for a = 

and in Proposition 1 of [47] for < a < 1) and the law, denoted M.a,ev, of M := ^XiVi 
is characterized by one of the following equalities: 

(i) for a = and 6* > 0, 

J Mo,eHdx){z-x)-'^ =cxp(^-e J i^{dx)\og{z~x)^ (zeC\R); (7.1) 

(ii) for ae (0,1) and 9>0, 

J Mc,.ev{dx){z - x)-'^ = (^j v{Ax){z ~ x)'^'^ ' (zeC\R); (7.2) 

(iii) for a G (0, 1) and 9 G (-a, 0), Ma,ev G V-e and (7.2) holds true; 

(iv) for a G (0, 1) and ^ = 0, MafiV G Vq and 

( Mafiv{<lx)\og{z~x)=a-^\og I v{dx){z-x)°' (zGC\R). (7.3) 

Thus the correspondence v i— )> Ma,ev defines a map from Va to V , for which we also write 
■M.a,e, by a slight abuse of notation. Since we are not aware of any reference in which 
intcgrability of the transformed measure Aia.ev claimed in (iii) or (iv) has been shown, 
the proof shall be given below by applying our result. As will be seen later (Proposition 
7.2(i)), such a property is also needed for further discussion. In the subsequent argument, 
we often use the following three equalities: 

/•oo 

/ dss'^-^e^'^e-"-^ =r(6l)(l + u)-^ (61 > 0), (7.4) 

/>oo 

Ca dss~°'-^e-''{e-'"' -l) = l-{l + uY (0<a<l), (7.5) 
Jo 

/■oo 

/ dss-ie-"(e^"'*-l) = -log(l + ?i), (7.6) 



JO 



where w G C is such that Rcu > — 1. 

Lemma 7.1. Lei a G (0, 1) and 6* G (— a, 0] . Then Ma,ev ^V-e for any v ^Va- 
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Proof. Define M* := ^ |Xi|T^. By the assumed independence, 



OO „ 

TT / Kd2;)e-^^'l-l 



is>0). (7.7) 



We will apply Theorem 3.2 by taking g{s) = J^i'(dx)eyi'p{—s\x\). For this purpose, with 
the help of Fubini's theorem, observe from (7.5) that for < a < 1, 

C f°° e"-^^ f 

dz^|g(z)-lH / Kdx)(l + A-i|x|)"-l<oo, (7.8) 

which implies that A* (5) = 0. Therefore, (3.8) holds for all A > and reads 
A'' 



1 



-e/a 



(7.9) 



i(^^Kdx)(l+A-Vir) 0<a<l, 0^0, 

--log / u(dx)(l + \-^\x\)°', O<a<l,0 = O. 

a Jr 

The desired intcgrability can be seen as follows. For 6 = 0, the left-hand side of (7.9) 
with A = 1 becomes — i?[log(l + M*)] by (7.6), while the right-hand side is finite by 
V G Va- This proves that Aia,Qi^ G Vq whenever v G Va- Similarly, for 6 G (— a,0), by 
using (7.5) with —9 in place of a, we see that the left-hand side of (7.9) with A = 1 
equals (£'[(1 + AI*)-'^] - l)/0. Therefore, Ma,ety e V-e is implied by ly^Va- □ 

Remarks, (i) In view of (7.9), the reader will be able to see that the dichotomy result 
of Feigin-Tweedie type [14] (originally shown for PD(0,6')'s) holds true for PD{a,9) 
with 6 > 0. More precisely, if v gV \ Va and M* is as above, then M* = 00 a.s. 
Thus the maximal domain of Ma,e with > is identified with Va- In the Dirich- 
let process case, Cifarelli and Regazzini [6] gave a proof of this fact by using (7.1). 
Their cutoff argument also applies for < a < 1, as seen in the following. Considering 
M* := ^ |Xi|l{|Xi|<n}l^i for each n = 1, 2, . . . , which corresponds to the image measure 
Vn of V under x 1— > 1 [_«,„] (cc), we have (7.9) with M* and f„ in place of Af* and v, re- 
spectively. Under the assumption that v ^ Va, letting n — > 00 yields E[{1 + AI* /X)^^] = 
and hence M* = 00 a.s., as long as > 0. 

(ii) The above proof makes it almost obvious that the Markov-Krein identity itself 
can be recovered from Theorem 3.2 with g(s) = 'ijj^{±s), where V'l' is the characteristic 
function of 1^. These choices are allowed for any v ^Va because for all x € R. 



00 



j.. e " \^^xz -^1 ^ / Aolog(l + |a;|)H-Bomin{|x|,l}, a = 0. 



"+!' ' - +S„min{|x|,l}, 0<a<l, 

where Aa and Ba are constants depending only on a. (For the proof of this inequality, 
first observe that the left-hand side is bounded above by |a;|r(l — a) and then show it 
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for x > 1 by, for example, dividing the doraain of integration into (0, l/x) and [1/.t,oo).) 
In addition, (7.5) and (7.6) show that A* (V',y(±-)) = 0- Therefore, analogously to (7.7), 



i^M^A^s) = £;[exp(±V^sA/)] = 
and we obtain, by (3.8), 



i=l 



(7.10) 



X' dss'-^e-^'{i'M^,,{±s)-l) = R^.ei^\-'' dz -^(V'.(±z) - 1) j (7.11) 

for aU A > 0. In the same way as before, (7.11) can be converted to (7.1)-(7.3) with 
z = ±\/—lX by making use of (7.4)"(7.6). 

The rest of the paper is devoted to describing various consequences of the above cal- 
culations and the Markov-Krein identity. 

Proposition 7.2. (i) Let Q < f3 < a < 1 and 9 > —(3. Then, as a map defined on Va, 

Mp,eoMa,-p=Ma,e- (7.12) 
(ii) Let < Q: < 1 and 9 > —a. Suppose that v G Va- Then, for a.e. t € R, 



e+an-l 



n=l ' "''^"2=1 "i \ j=l I 



(iii) Let < a < 1 and 9 > —a. Then Ma.e defined on Va is injective. 

Proof, (i) First, note that by Lemma 7.1, the composition Aip^e o Aia,-i3 is weh defined 
on Va under the conditions assumed. (7.12) with 9 = becomes Mpfi o Ma.-i3 = A^a,o, 
the proof of which requires careful handling of branches. Taking an arbitrary v G Va and 
putting = {AijB.o o Aia.-p)'^, we observe from (7.2) and (7.3) that 



R 



i^'(da;) log(z — x) = — log / Ma.-i3v{dx){z — x)^ 



loe 



R 



13/a 



i'{dx){z — x)"' ] = — log exp — log / v{dx){z — x) 



R 



It should here be noted that since < /3/a < 1, 



arg exp 



-1— arg / i'{dx){z — xY 
" Jr 



arg / v{dx){z — x)". 



R 
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Combining these equalities yields 

/ z^'(dx) logfz — a;) = 7- • — log / h'(dx)(z — x)" ^ [ Ma o'^idx)\og(z ~ x), 
Jr P a Jr ' 

where the last equality follows from (7.3). Since z € C \ R is arbitrary, this implies 
v' = M-afiV, as required. It only remains to deal with the case where < /3 < a < 1 and 
6* 7^ is such that > — /?. Calculations are quite similar to those in the previous case 
and so are left to the reader. The proof of (i) has already been done. 

(ii) By virtue of (7.11), (7.13) is a direct consequence of Lemma 3.1(i). 

(iii) First, suppose that MafiV = Mafiv' for some £ Va- Then, by (7.3), 
^^v{dx){l + v^A-ix)" = /j^z/'(dx)(l + V~A-ia;)" for all A G R\ {0}. This, together 
with (7.5), implies that v ~ v' . Next, let 0^0 be such that d > —a and assume that 
■A^Q.gi^ = M-afiv' for some v' G Va- Taking care with branches, we see by (7.2) that for 
each A e R \ {0}, there exists an integer n(A) such that 

log/ v{<^x)(\^^J~^\-^xf -\o%\ v'{Ax)(\^^f^\-^x)"- ^2'Kn{\)ad-^. (7.14) 
Jr Jr 

However, since the left-hand side is continuous in A and tends to as |A| — ?> 00, we have 
n(A) = 0. This proves v = v' and therefore Ma.d is injective. □ 

Let us make some comments on Proposition 7.2. The assertion (i) extends Theorem 
2.1 of [25], where (7.12) with /3 = is shown and applied for a sampling procedure. We 
can also understand (7.12) via (7.11) combined with R^^ ° Ra,-i3 = Ra,e- In general, a 
random variable M having {Ma,i3 ° -M-y^s)!^ distribution (if such a law is well defined) is 
constructed by M = J2jiJ2i-^ij^ij''^^Wj°''^\ where {Xy} is a family of i.i.d. random 
variables with common law I/, (T^/"''^^)°li has PD(a,/3) distribution and {v}^-^^)°Zi, 

{^i2'^^)'i^i^--- ^^"^ PD(7,(5)-distributed. Moreover, these random elements are required 
to be mutually independent, (ii) and (iii) are motivated by some results of Lijoi and 
Regazzini [36] concerned with A4o,e- Applying to (7.11) a general complex inversion for- 
mula for Laplace transforms (the Bromwich integral), the authors provided an expression 
for the characteristic function of Aio^giy. It is not clear how to verify equivalence between 
such a formula and (7.13) directly. They also proved a prototype of (iii), that is, that 
■Mo^e '-Vq^V is injective for each 6* > 0. 

As observed [26] (Remark 4.2) in a much more general scheme called a species sampling 
model, Cauchy distributions arc fixed points of every Aia.e- Alternatively, this also follows 
from (7.10). It should be emphasized that any property of PD(q;,0) as a law on Voo is 
irrelevant to the validity of this fact, except the one which ensures that ^ < 00 

a.s. The study of the inverse problem would require a deeper understanding of Aia,e 
itself or PD(a,0). In this respect, Yamato [55] gave a partial result concerning the case 
a = 0. It was shown in [36] that for any v (IVq, Mq^iu = v occurs only when v is Cauchy 
or degenerate. It seems reasonable to expect that such an equivalence would hold true for 
any Ma, 6 and we could approach this problem by the use of results presented in Section 
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3. But here we shall show only a minor extension by mimicking the proof of Theorem 5 
in [36], that is, by reducing to an ordinary differential equation via the Markov-Krein 
identity. 

Proposition 7.3. Let < a < 1 and v e Va- Then M-aS-aM ~v ij and only if v is 
Cauchy or degenerate. 

Proof. Supposing that Ma^i-aV = v, we only have to show that v is Cauchy or degen- 
erate. By combining this equality with (7.2), we easily see that f±{t) := ^^v{dx){t ± 
y^x)" solves the equation a-^f±'{t) = /±(i)^~^/"(t > 0) and therefore f±{t)^^°' = 
t-|-/±(0)^/"(t > 0). Here, since /+(0) and /-(O) arc complex conjugate with | arg/±(0)| < 
7ta/2, we have | arg(/±(0)^/")| < 7t/2 so that /±(0)^/" = cr^ ± ^f^m for some cr^ > q 
and m S R. Consequently, f±{ty^" = t + ± \/—lm{t > 0). Furthermore, by noting 
the continuity of arg/±(t), we get f±{t) ~ {t + ± ^/^m)" for every t > 0. With 
the help of (7.5), we conclude that v is a Cauchy or degenerate distribution with 
■ip^{t)^cx.p{y/^mt-a'^\t\). □ 
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